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CHAPTER  1 


Prologue 


Study  Objective  and  Methodoloc 


The  primary  objective  of  this  investigation  was  to  analytically 
determine  tne  effect  of  a  plane  dielectric  interface  upon  the  radiation 
loss  of  a  curved  channel  dielectric  waveguide  embedded  in  a  dense  substrate 
and  bent  in  a  plane  parallel  to  that  of  the  material  interface.  The  radiation 
loss  from  bent  dielectric  waveguides  is  of  interest  for  a  variety  of  potential 
applications  in  integrated  optics.  Knowledge  of  the  factors  influencing  the 
radiation  from  bend  dielectric  waveguides  will  allow  the  design  of  integrated 
optical  component  characteristics  which,  in  one  application,  may  result  in 
decreased  radiation  to  conserve  signal  energy  or  which,  in  another  applica¬ 
tion,  may  achieve  enhanced  radiation  effects  in  order  to  optimize  the  coupling 
between  vatious  optical  waveguides.  An  obvious  example  for  which  the  conser¬ 
vation  of  signal  energy  is  of  prime  concern  is  the  dielectric  fiber  trans¬ 
mission  line  used  in  simple  point-to-point  data  transfer  applications. 

Examples  of  integrated  optical  components  in  which  some  specific  level  of 
radiation  is  desired  include  optical  filters  (e.g.,  ring  resonator  discrimi¬ 
nators;  cf.  [1])  and  optical  directional  couplers  (cf.  [2]).  Figure  1.1 
illustrates  potential  realizations  of  these  latter  components. 

Arnaud  (cf.  [3])  has  treated  radiation  loss  from  a  slab-loaded,  bent 
dielectric  rod,  and  this  structure  is  the  one  described  in  the  literature 
that  has  the  most  in  common  with  the  configuration  considered  herein. 


(a)  A  RING  RESONATOR  DISCRIMINATOR 


(b)  A  DIRECTIONAL  COUPLER 


Figure  1.1  Dielectric  Channel  Waveguides  in  Potential 
Integrated  Optics  Applications 
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However,  primarily  because  of  the  finite  extent  of  the  slab  and  the  fact 
that  the  refractive  indices  of  the  slab  and  the  rod  are  equal,  Arnaud's 
problem  differs  significantly  from  that  discussed  in  this  report.  Indeed, 
the  differences  are  substantial  enough  to  preclude  extensive  corroboration 
ot  either  study's  results  by  direct  comparison. 

Figure  1.2  illustrates  the  fundamental  geometry  of  the  problem  and 

depicts  the  coordinate  reference  frames  chosen  for  the  analysis.  The 

refractive  indices  of  the  channel  waveguide,  the  substrate  material,  and 

the  covering  material  nave  the  values  n,  n^,  and  n^  ,  respectively. 

These  refractive  indices  are  ordered  according  to  the  relation  n  >  n^»  n.  , 

the  relation  n2»n^  expressing  what  is  meant  by  "dense"  as  applied  to 

the  substrate  material.  A  global,  circular-cylindrical  reference  frame 

with  spatial  coordinates  r,  $,  and  z  has  its  origin  at  the  point  0,  and 

a  local  cartesian  reference  frame  with  spatial  coordinates  x,  y,  and  s 

has  its  origin  at  the  point  (Rq,  0,  0)  relative  to  the  global  frame.  The 

extent  of  the  channel  waveguide  is  defined  by  the  radial  boundaries  at  r" 

and  r"1  and  by  the  coordinates  z  =  0  and  z  =  z",  which  define  horizontal 

planes  between  which  the  guide  is  confined.  The  quantity  Ro  is  given  by  the 

expression  R  =(r'+r'"]/2  and  is  defined  to  be  the  radius  of  curvature  of  the  channel 
waveguide. 

A  propagating  mode  with  phase  factor  exp  i(wt  -  kQy  RQ<?)  is  assumed 
to  exist,  where  kQ  =  w/pQc q  is  the  free  space  propagation  constant  and  y 
is  the  normalized  (to  kQ),  complex  propagation  constant  of  the  mode  in  the 
curved  channgel  waveguide.  Given  that  the  normalized  propagation  constant 
for  the  straight  channel  waveguide  is  y  ,  the  analysis  proceeds  to  compute 
a  first  order  representation  of  y  of  the  form 

y  =  y  -  ia  ,  a  >  0 
o 


(l.D 
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wnere  oc  Is  the  normalized  attenuation  coefficient  of  the  mode.  The  moce 

attenuation  due  to  continuous  radiation  from  the  bend  following  propagation 

through  an  angular  sector  <?  of  the  channel  guide  is  then  given  by  the 

factor  exp  (-ak  R  p) . 

oo 

Once  a  has  been  determined,  the  study  objective  is  attained  via  a 
comparison  of  a  with  the  normalized  attenuation  constant  of  a  bent  channel 
waveguide  embedded  in  a  homogeneous  medium. 


where  >  0  and  Q(0)  <  0  are  parameters  depending  upon  material 
characteristics  of  the  guide  (including  physical  dimensions)  and  its 
surrounding  medium  as  well  as  the  specific  mode  assumed  to  be  propagating 
in  the  guide.  As  a  result  of  this  investigation,  the  normalized  atten¬ 
uation  constant  a  for  the  bent,  rectangular  dielectric  waveguide  embedded 
in  a  dense  substrate  below  a  material  interface  is  found  to  be  of  the  form 


where  C.  >  0  and  q^(0)  <  0  depend  upon  the  channel  guide,  substrate, 
and  covering  medium  characteristics  (including  the  physical  dimensions  of 


the  guide)  and  upon  the  specific  mode  assumed  to  be  propagating  in  the 


guide. 


Three  observations  are  of  importance.  First,  the  quantities 

and  Ci  are  in  general  unequal,  as  are  Q(0)  and  q2 ( 0)  (although  a 

certain  formal  similarity  between  the  latter  two  exists).  Second,  the 

fundamental  exponential  dependence  of  both  a'  and  a  on  kQ  is  tne 

same.  Third,  the  algebraic  dependence  of  a*  and  a  on  R  is  distinctly 
i  3/2 

different.  The  (5—)  algebraic  dependence  of  a  upon  R  arises  as 
Ko  0 

the  leading  term  of  the  asymptotic  expansion  of  a  and  thus  imparts  a 

lateral  wave  characteristic  (cf.  4,  Section  5.5)  to  the  radiation  loss 

of  a  rectangular  dielectric  guide  in  a  aense  substrate  below  a  material 

interface.  Indeed,  the  absence  of  a  („ - )  leading  term  in  the  asymDtotic 

0 

exransion  of  a  may  be  thought  of  as  resulting  from  the  first  order  can¬ 
cellation  of  source  and  image  field  due  to  the  presence  of  the  interface. 

C_ 


As  a  means  of  comparing  the  influence  of  the  factors 


and 


- 1— : y.-y  on  the  radiation  loss  of  their  associated  waveguiding  structures, 

w ' 

a 1  and  a  are  plotted  versus  kQRo  for  bent,  rectangular  dielectric  wave¬ 
guides  whose  dimensions  and/or  refractive  indices  are  chosen  so  that 
Q ( 0 )  =  q^ ( 0) .  This  not  only  assures  that  the  exponential  variation  with 
Rq  will  be  the  same  in  each  case,  but  it  also  corresponds  to  a  situation 
in  which  propagation  constant  along  the  central  axis  of  each  guide  is  the 
same.  For  the  cases  considered,  a  is  found  to  be  always  less  than  a' 
and  is,  in  particular,  less  than  an  order  of  magnitude  smaller  than  a' 

for  k  R  =  500  and  between  one  and  two  orders  of  magnitude  smaller  than 
0  0 

a'  for  k  Ra  =  5000.  This  indicates  that  one  may  in  some  instances 
0  0 

achieve  a  significant  reduction  in  the  radiation  loss  of  a  bent,  rectangular 
dielectric  wavequide  by  embedding  it  below  a  high  contrast  dielectric  inter¬ 
face  in  the  material  of  greater  refractive  index. 


7 


1 . 3  Analytical  Preliminaries 

Inis  report  will  not  present  a  derivation  of  the  general  form  of  a, 
but  will  merely  quote  the  result  as  given  in  [5],  Chapter  9.  Thus 


where 


p  =  2  j  VEo  *VS 

oo 

C  =  [  {£"  x  6H  -  6E  x  H-}  •  a  d 

Jo  o  y  z; 


/  i  a  V 

\  I  -2; 


(1.3) 


0.4) 


In  (1.3),  the  integration  extends  over  an  infinite  plane  containing  a 

cross  section  of  the  channel  waveguide,  and  yQ  is  some  point  chosen 

outside  the  bend  (beyond  r‘")  of  the  curved  channel  guide  where  the  fields 

have  decayed  sufficiently  (cf.[2j).  The  fields  Eq  and  Hq  are  those  of 

the  straight  channel  guide,  and  E~  and  H~  are  the  fields  of  the 

straight  channel  traveling  in  the  "negative''^  (or  s)-direction.  The  fields 

6E  and  6H  constitute  that  portion  of  the  fields  outside  the  curved 

channel  guide  which  have  been  reflected  from  the  caustic  or  turning  point 

beyond  which  the  mode  must  radiate  in  a  radially  outward  direction.  The 

fundamental  analytical  task  of  this  investigation  was  that  of  expressing 

tne  fields  TE  and  6H  in  terms  of  known  quantities,  i.e.,  in  terns  of 

E  ,  H  ,  y  ,  R  ,  etc. 
oooo 

As  a  result  of  the  circular-cylindrical  symmetry  of  the  global 
reference  frame,  our  subsequent  analysis  will  be  facilitated  by  concentra¬ 
ting  on  the  field  components  Ez  and  hz,  although  we  must  eventually 
consider  all  of  the  field  components.  Via  Maxwell's  equations  and  standard 


mwmmmrnmmmmimm 


analysis  for  circular-cylindrical  geometry,  we  find  that  we  may  express  tne 


four  remaining  field  components  in  terms  of  Ez  and  whence 
,2  3H  R  3E 


<fc2 '7'V  «oko  77  -  lY  k-  ^  :'Z 


3z 


or  3z 


o  *2  o  R 

(IT  +  -*-*)  E  =  -  (  y  k  -  H  + 

'  'f  r\  '  r\  Y'  7 

3  z 


.  Y  N  -  i  1  1  ^  U 

o  o  r  z  or3z  z 


0  ~,2  .  2  3E  R  ~u 

<k  +  &  "♦  =  -TT  ko  77  -  iYko  7*  77 


,  ~2  2  9  R  a2 

+  i--  lu  =  !L_  Yk^  —  E  +  — —  H 

U  a,2  J  r  C  Y  o  r  z  3r3z  z 

3z  o 


(1.5a) 

(1.5d) 

0.5c) 

(1 .5d) 


- -  2  2  2  2  2  2 

where  ;  =  At  /e  ,  k  =  n,k  for  z  >  0,  k  =  n-k  for  z  <  0. 

ooo  1  o  2  o 


We  will  show  in  Chapter  2  that  both  E^  and  Hz  have  Fourier-Bessel 


representations  of  the  form 


E_(r,<{>,z)  =  e'1v$ 


Ez(a,z)Jv(kQar)adot 


(1.6a) 


Hz(r,<?,z)  =  e 


=  q-‘M> 


Hz(a,z) Jv(kQar)ada 


(1.6b) 


All  of  the  desired  field  components  have  similar  representations,  and 
the  z-variation  is  such  as  to  yeild  the  following  equivalence  of  operators 


<-**  *  S)  ■  A2 

3z  0 


0.7) 


for  the  field  transforms. 


For  example,  consider  the  electric  field  component  E^  given  by 


S  =  e'"*  I  E^(a.2)Jy(fc0«r)adc 


(1.8) 


.... - - 


1 
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g u  2  - 

E,(a,z)  is  the  field  transform  of  E^  and  ( — ^  +  k  )E*(a,z)  = 
2  9-  V  32 

kVE.. 

o  4> 

Substitution  of  (1.6a),  (1.6b),  and  (118)  into  (1.5a)  and  making  use  of 
(1 .7)  gives 


i;  a 

o 


J ' (k  ar) 
vv  0  ' 

J  (k  a  )* 
vv  o  r' 


—  E 

32  2 


(1.5a) 


In  a  similar  fashion  we  find  the  following  as  well: 


a2E 


c  yR  _  J  ‘  (k  ar)  . 

o  o  u  +  a_  vk  o  3_  f 

r  2  k  J  (k  ar)  32  2 

o  v  o 


2~  JL  u  i  *  nl,2  Jv  (koar)  F 

a  0  "  '  rk  32  H2  "  r  Jjk  a'r)  z 


o 

2 


v'  o 


2  p,  .  yRonl,2  p  a  Jv(koar^  3  fj 

a  ^  ~W ^  »  z 


(1.9b) 


(1.9c) 


(1.9d) 


The  relationships  (1.10a)  through  (l.lOd)  between  the  field  (Fourier-Bessel ) 
transforms  will  be  of  use  in  Chapter  2. 
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CHAPTER  2 

INTEGRAL  REPRESENTATIONS  OF  THE 
CARTESIAN  FIELD  COMPONENTS 


2.1  Introduction 

This  chapter  is  devoted  to  finding  integral  representations  for  tne 
Cartesian  field  components  E^  and  Hz-  As  mentioned  in  Chapter  1, 
knowledge  of  these  allows  one  to  determine  the  remaining  four  field  compo¬ 
nents  via  (1.5a)  through  (1.5d).  The  components  Ez  and  Hz  are  pre¬ 
ferred  by  the  circular-cylindrical  geometry  of  our  problem  to  the  extent 
that  they  both  satisfy  the  homogeneous  (in  source  free  regions),  scalar  wave 
equation,  i.e.,  the  Helmholtz  equation.  This  fact,  together  with  the 
assumed  form  of  the  field  variation  with  the  spatial  coordinate  $> ,  allows 
specification  of  the  Fourier-Bessel  transforms  of  Ez  and  Hz  to  within  mul¬ 
tiplicative  factors  which  must  be  determined  through  the  application  of 
boundary  and  source  conditions.  Explicit  evaluation  of  these  is  the 
subject  of  Section  2.3. 

Although  at  first  glance  the  final  form  of  the  integral  representations 
of  Ez  and  Hz  seems  too  cumbersome  for  practical  use  due  to  the  complexity 
of  the  required  integrations,  it  will  be  seen  in  Chapter  3  that  the  integral 
representations  derived  in  this  chapter  are  amenable  to  asymptotic  analysis 
when  kQR0  »  1,  i.e.,  when  the  radius  of  curvature  of  the  bend  channel  guide 

is  large  relative  to  the  wavelength  of  the  propagating  fields. 
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2 . 2  Extended  Fourier-Bessel  Representations  for  Field  Solutions  of  the 
Helmholtz  Equation 

Let  F(r,q,z)  =  4>(r ,z)e’lv^  be  a  scalar  field  satisfying  the 
Helmholtz  equation  for  source-free  regions, 

2  ? 

(V  +  k  )F(r,<?,z)  =  0,  then  for  v  =  k0YRQ> 


o  2.2d2 

,ijl  r|_  + j>5  +  k2  - 

r  3r  3r  ^2  r2 


$(r,z)  =  0 


(2.1) 


2,2 


2  .  J.,Z 


where,  as  before,  k  =  n-jkQ  for  z  >  0,  and  k  =  n^k^  for  z  <  0. 

In  this  section,  o(r,z)  is  a  generic  symbol  representative  of  eitner 
Ez(r,z)  or  Hz(r,z),  since  the  mathematical  manipulations  required  herein 
do  not  require  us  to  distinguish  these  two  field  components.  Appropriate 
associations  between  Ez(r,z)  and  H  (r,z)  and  the  results  of  this  section 
will  be  made  as  required  in  later  sections. 

From  the  form  of  F  (r,$,z)  we  know  that  there  is  a  Fourier-Bessel 
transform  representation  (cf.  [6])  for  $(r,z),  whence 


$(r,z)  =  |  i  (a,z)J  (k  ar)ada 


v  o 


(2.2) 


o 

00 

<5>(a,z)  =  k^  jt>(r, z)Jv(kQar)rdr  (2.3) 

o 

where  is  a  Bessel  function  of  order  v.  We  shall  refer  to  <t>(r,z) 
as  the  field  and  to  <£(r,z)  as  the  field  transform. 

Substitution  of  (2.2)  into  (2.1)  produces  the  result 


3  z 


2  *  k|(nl2 


%(a,z)  *  0 


(2.4) 
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where  n‘  =  n-j  for  z  >  0,  and  n'  =  x\p  for  z 
If  we  consider  z  =  z',  where  z"  <  z‘  <  0, 
sources  lie,  and  if  we  likewise  require  that  the 
S  zj  -*00,  then  (2.4)  has  solutions 

-k  u.z 

$  (a, z)  =  A(ot)e  0  ,  z  _>  0 


<  0. 

as  a  plane  in  which  point 
field  should  vanish  for 


(2.5) 


-k  u9z  k  u.z 

4>v  (a,z)  =  Bp(a)e  0  +  Bn(a)e  z‘  <  z  <_  0 

k  u  z 

5V  (a,z)  =  C(a)e  ,  z  <  z* 


(2.6) 

(2.7) 


where 

u1  =  (a2  -  n2)'2  ,  Re  u]  >  0  (2.8) 

Up  =  (a2  -  r\p)^  ,  Re  Up  >  0  (2.9) 


With  (2.5),  (2.6),  and  (2.7)  we  have  determined  the  field  transforms 
except  for  a  set  of  multiplicative  factors  dependent  upon  a.  In  the 
following  sections  we  will  apply  source  and  boundary  conditions  to  deter¬ 
mine  these  factors.  At  that  time  it  will  be  necessary  to  identify  our 
field  4>(r,z)  as  either  the  electric  field  Ez(r,z)  or  the  magnetic 
field  Hz(r,z). 

To  facilitate  the  asymptotic  analysis  of  Chapter  3,  it  is  desirable 
to  extend  the  range  of  integration  in  (2.2)  to  the  negative  a  axis.  To 
do  this  we  recall  that 

Jv(oa)  =  Ju(H^(oa)  +  H^(aa)} 

and  employ  the  circuital  relations  for  the  Bessel  functions  so  that  (2.2) 


becomes 


«(r,z)  =  h 


OO0 


-ITT 


0  (a,z)F^(k  ra)cxda 
v'  v  o 


(2.10) 


Figure  2.1  illustrates  the  integration  contour  of  0(r,z)  in  tne  complex 
a-plare.  The  branch  cuts  associated  witn  n,  and  n„  are  those  which  must 

i  2 

not  be  crossed  if  we  are  to  insure  that  Re  u,  >  0  and  Re  u„  >  G.  Some 

i  L 

slight  loss  is  assumed  in  both  the  covering  medium  and  the  substrate  in 
order  to  impart  a  small  imaginary  component  to  both  n^  and  r^.  The  branch 
cut  running  along  the  negative  real  axis  insures  that  %  "(k^ar)  in  (2.10) 
is  single-valued. 


2 . 3  Evaluation  of  the  Field  Transforms 

2.3.1  Equivalent  Polarization  Sources 

The  concept  upon  which  the  analysis  of  this  section  rests  is  that  of 
equivalent  polarization  currents  (cf.  [7])  which  can  be  used  to  convert 
the  sourceless  three  medium  problem  into  an  equivalent  two  medium  problem 
with  sources  below  the  dielectric  interface.  This  approach  attributes  the 
fields  existing  outside  the  dielectric  channel  to  the  radiation  from  an 
array  of .polarization  current  sources  contained  within  the  geometrical 
boundaries  of  the  channel  waveguide  but  embedded  in  the  substrate  medium 
only.  This  circumstance  is  illustrated  in  Figure  2.2. 

The  justification  for  this  procedure  follows  from  a  consideration  of 
the  differential  form  of  Ampere's  Law,  i.e.,  the  curl  H  relation  from 
Maxwell's  equations.  In  the  channel  region,  we  have 

V  x  h  =  i w  n (2.11) 


and  in  the  substrate  region  the  fields  must  satisfy 


(a)  THE  ORIGINAL  THREE  MEDIUM  PROBLEM 
WITHOUT  SOURCES 


(b)  THE  EQUIVALENT  TWO  MEDIUM  PROBLEM 
WITH  POLARIZATION  CURRENT  SOURCES 


Figure  2.2  Alternative,  Equivalent  Problem  Representations 


16 


V  x  H  =  i  w  n|e  E  (2.12) 

2  _ 

In  (2.11)  we  add  and  subtract  jwn2£QE,  then  (2.11)  becomes 

V  x  H  =  i  w  (n2  -  n2)  cqE  +  i  w  n2eQE  (2.13) 

We  now  identify  the  first-term  on  the  RHS  of  (2.13)  as  the  equivalent 
polarization  current  density  Jp,  thus 

Jp  =  i  w(n2-  n2)  i:o£  (2.14) 

and  (2.13)  is  now 

V  x  h  =  J +  iwn^E  (2.15) 

Comparison  of  (2.15)  with  (2.12)  shows  that  we  may  now  characterize 

the  fields  below  the  dielectric  interface  through  the  use  of  (2.15)  alone, 

as  long  as  we  keep  in  mind  that  the  equivalent  polarization  currents  exist 

only  within  the  region  occupied  by  the  channel  guide.  Since  the  only 

refractive  index  that  appears  explicitly  in  (2.15)  is  n£,  our  problem  is 

seen  to  be  equivalent  to  that  of  a  two-medium  interface  problem  with 

sources  as  illustrated  in  Figure  2.2b. 

It  should  be  emphasized  again  that  to  solve  the  problem  discussed 

herein  we  assume  that  the  fields  E  and  H  in  the  channel  guide  are 

oo 

known,  and  the  implication  of  this  assumption  in  the  present  discussion  is 
that  Jp  will  be  given  explicitly  as 

Op  =  i  w(n2  -  n2)coEo  (2.16) 

We  have  thus  converted  information  about  known  channel  fields  into  a  known 
distribution  of  equivalent  current  sources.  Since  Eq  will  in  general  be 
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capable  of  resolution  along  each  of  the  three  mutually  perpendicular  axes 
of  the  local  coordinate  system  (cf.  Figure  1.2)  (2.16)  reveals  that  we 
will  in  general  have  polarization  currents  directed  along  all  three  of  the 
local  coordinate  axes.  The  remainder  of  this  chapter  is  devotee  to  finding 
the  field  transform  amplitudes  associated  with  each  unique  component  of  the 
polarization  currents.  The  total  electric  and  magnetic  fields  outside  tne 
channel  will  ultimately  be  constructed  via  superposition  of  the  results 
obtained  for  the  special  cases  to  be  considered  currently. 


2-3.2  z-Directed  Phased  Arrays  of  Polarization  Currents 
We  consider  a  z-directed,  phased  array  of  polarization  currents 
located  at  a  radius  r1  from  the  z-axis  and  a  distance  z'  below  the 
dielectric  interface  as  illustrated  in  Figure  2.3.  Analytically,  polar;- 
zation  line  current  densities  of  unit  amplitude  along  the  z-direction  are 
given  by 

Jp  =  az 6  (z-z‘ )  & (r-r' )e-i (2.17) 

hence  reference  to  (2.16)  gives 

Jp  =  iw(n2-n2)eoEoz3pdr'dz‘  (2.18) 


To  facilitate  the  analysis  of  this  section  we  shall  consider  a 
z-directed  electric  type  Hertzian  potential  of  the  form 


n  e(r,<j>,z)  =  a2e'lv<J> 


oo 

n 2(a,z)0v(kora)ada, 


so  that 


n®(r,z) 


n 2(a,z)Jv(kQra)ada 


(2.19a) 


'  - . . - 
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and 

oo 

n®(a,z)  =  k2  J  H®(r,z)Jv(kQra)rdr  (2.19b) 

o 

In  this  section  we  will  be  dealing  with  the  cartesian  component  the 
electric  field  only,  because  H  =  0  for  a  z-directed  Hertzian  potential. 

The  Hertzian  potential  satisfies  the  inhomogeneous  wave  equation 


(cf.  [8]) 


(V2  +  k2j2)ne  = 


-J 


iweon2 


(2.20) 


where 


kl  2  =  nl  2ko’  thuS  us^ng  (2.17)  in  (2.20) 


(V2  +  k2  «)n ®(r,z)  =  -  ( - - — *-)  <S  (z-z1)  6  (r-r*) 

1  \  UI  c  r\C 


(2.21) 


i  w  e  n_ 
o  2 


Substitution  of  (2.19a)  into  (2.21)  and  use  of  the  Fourier-Bessel 
representation  of  6  (r-r')  yields 

y  -  &?.  J  *;<«.*)  -  [-^]  vk„r,“> 6  (z-z,)  (2-22) 

*•  *  o  2 

2  2  ^ 

where,  as  in  (2.8)  and  (2.9),  u^  2  =  (a  -  2) 

Following  the  same  line  of  reasoning  as  in  Section  2.2,  we  have  the 
following  solutions  for  n®(°i,z)  from  (2.22): 

-k  u.  z 

fi  ®(a,z)  =  A(a)e  01  ,  z>0  (2.23a) 

-k  u~z  k  U-Z 

S®(a,z)  =  Bp(a)e  +  Bn(a)e  0  L  ,  z’<  z<0  (2.23b) 

-k  u_z 

n®(a,z)  =  C(a)e  0  1  ,  z  <  z’  (2.23c) 
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We  solve  for  the  transform  coefficients  A(a),  Bp(a),  Bn(a),  and 
C(ot)  via  application  of  source  conditions  at  z  =  z‘  (two  equations)  and 
boundary  conditions  at  z  =  0  (two  equations);  this  yields 


u2 V‘  fW"‘>?  A>U2Z' 

2  2  2  (  e 
uln2  +  u2nl  ^-1u2W£on2  J 


(2.24a) 


k  J  , ( k  r‘*)r‘  k  u.z' 

6  JO  ■  - e  “  2 

p  i2u0we  n0 

2  0  2 


(2.24b) 


„  ,  ,  ko')v<kor'a)r'  f“2nl-ul"2  1  .koV' 

B  (a)  =  - 2 -  1  ~'2 - 2f  e 

i2u2wcQn2  (u2n1+n1u2j 


(2.24c) 


k  J.(k  r'a)'.-'  f  -ku,z'  u0n^-u,n,  kuz'i 

CM-  -°-v— < — e  02  +  1  I  -  e  0  2  1  (2.24c) 

i2u2wEon2  l  Wuln2  J 


u2n.j  +  u^n2 


We  can  now  use  (2.24a)  through  (2.24d)  in  (2.23a)  through  (2.23c)  in 
order  to  produce  the  following  expressions  for  the  potential  transforms 


=  e,  .  Vv(kor'a)r'  f  2V 2 

ji  (oi,z)  2  <  2  2 

i2u2W£Qn2  (u1n2  +  u2n1 


ko(V'-“lz) 


z  >  0 


(2.25a) 


n  |(cx,z) 


ko'Mk0r,a)r'  [  -koU2 ' Z_Z ’  u2nl  ”  Uln2  koU2^Z+Z,) 
2  <e  +  2  '2  e 

i 2u2  w  eQn2  (  u^  +  u^n2 


z  <  0 


(2.25b) 
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Substitution  of  (2.25a)  and  (2.25b)  into  (2.19a)  produces  the  desirec 


integral  representations  for  the  Hertzian  potential  in  the  covering  region 


(z  >  0)  and  the  substrate  region  (z  <  0): 


32tr,z)  -  -H  r'  (— - -2-)J  (k  r'a)J  (k  r-ale 

u,n0  +  u„n. 


k  (u0z'-u,z)  , 

o  2  i  a<xi 


o  12  2  1 


z  >  0 


(2.26a) 


VjI2-2!  ru?nruin2i  knu2(2+2'>i 


2nz  i  1  Lu  nr  +  u,n);- 


2  1  12 


0  (k  r‘a)J  (k  ra)~  ,  z  <  0  (2.26b! 


v  o  v  o  u 


where  we  have  used  ^  ,  the  intrinsic  wave  impedance  of  free-space. 

Proceeding  as  was  done  to  extend  the  range  of  integration  of  (2.2)  to 
yield  (2.10),  we  may  extend  the  range  of  integration  in  (2.26a)  and  (2.26b) 


to  give 


.  -ic  r  r 

nf(r,z)  =  -  -s9  -  f 


2  J 

°°e 


I - ~ H^(k  ar)J  (k  r'a)e  0 

-i  tt  u ^ n 2  +  u2n2  v  0  v  0 


ko(u2z'"ulz) 


z  >  0 


(2.27a) 


n  “(r,z) 


2  , 

4n0  J -i  tt 
2  °°e 


.2  .  2 


.  -V  !  f.‘kou2|z-2'l  .  u2nTuln2  .kou2(z+z,) 

“  a  (  6  >  r\  ri  S 


2  2 
u2ni  + u^n2 


H<2>(kom)0v(kor'a)  ,  Z£0 


(2.27b) 


4 


2.3.3  r-Directed  Phased  Arrays  of  Polarization  Currents 


We  consider  an  r-directed  phased  array  of  polarization  currents  located 
at  a  radius  r'  from  the  z-axis  and  at  a  distance  z'  below  the  dielectric 
interface  as  illustrated  in  Figure  2.4.  Figure  2.4b  is  a  cross-sectional 
view  while  Figure  2.4a  is  a  view  as  seen  from  the  covering  medium  looking 
down  on  the  z-axis. 


Analytically,  the  polarization  current  sources  of  unit  amplitude  are 


given  by 


hence 


5  =  a  6  (z-z1 )  6  (r-r' 
p  r 

T  •  iw(n2-^)eoE  3  dr'dz' 


(2.28) 


(2.29) 


The  approach  taken  in  this  section  is  to  consider  directly  the 
cartesian  components  of  the  electric  ana  magnetic  fields.  We  have  seen 
in  Section  2.2  that  and  Hz  satisfy  the  Helmholtz  equation  in  source- 
free  regions  anu  thus  that  their  associated  field  (Fourier-Bessel )  trans¬ 
forms  satisfy  equations  (2.5),  (2.6),  and  (2.7).  In  light  of  the  results 
in  Section  2.3.2,  we  will  hypothesize  from  the  start  the  following  forms 


for  the  field  transforms 


=  A®’m(a)e 


-ko(Uiz-u2z') 


z  >  0 


^’%,z)  =  B^m(a)e 


k  u_(z+z' ) 


*  B?-%)e  0  2 


-k  u.(z-z') 


z'  <  z  <  0 


(2.30a) 


(2.30b) 


k  U..Z 
0  2 


z  <  z' 


(2.30c) 


where  we  will  make  the  associations  A  §z  and  4™  =  Hz, 


rid  the  subscript 
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r  indicates  that  these  are  field  transforms  associated  with  r-directed 
polarization  currents. 

In  order  to  compute  the  eight,  as  yet  unknown  coefficients  of  the 
field  transforms  in  (2.30a)  through  (2.30c),  we  must  determine  anc  apply 

-v-  p  —  m 

the  boundary  and  source  conditions  applicable  to  4>r  and  .  This  process 
is  presented  in  detail  in  Appendix  A  and  yields: 


m  2u0 

Am  =  - Bm 

r  u-j  +  u2  p» 

„  u0  -  u,  „ 

gm  __  2 _ 1  gm 

nr  u^ +  u^  pr 

Ae  _  2u2n2 

r  "  2  2 

u2nl  +  Uln2 

2  2 
ge  _  u2nl  "  UIn2 

nr  u2nf  +  U,t>f 


u2  -  u1 

k  u„z 

-k  u„z 

o  2 

„  o2 

uT  +  u2 

e 

-  e 

2  2 
u^-j  -  u^ 

k  u„z' 

-k  u„z 

„  0 

n  o2 

2  2 
u^  +  u^2 

e 

-  e 

where  B  and  B*;  are  given  by 


iyR  k^ 

UK  «r') 


ar'5  K 

-tf"  ^(k°“r,) 


(2.31a) 

(2.31b) 

(2.31c) 

( 2 . 31  d ) 

( 2 . 3 1  e ) 

(2 .31 f ) 


(2.32a) 

(2.32b) 


Substitution  of  (2.31a)  through  (2.31f)  into  (2.30a)  through  (2.30c) 
yields  the  following  forms  for  the  field  transforms: 


fry*}™ 


5  ®(a,z) 


4  *(a,z) 


2u2n2  -k0(VV  }  .  .  n 

V  2~ - 2  e  *  2-° 

u2n1  +  u1n2 

Be  f  vi_LVLekou2(2+2,i 

Pr  l  u2nl  +  vl 

-k  u  j z-z  ‘  1  7 

+  sgn[(z-z' )]e  02  >  z  <  0 


(2.33a) 


-  m,  \  jvi 

4  r(a,z)  =  B 


4  J!(a,z)  =  B 


2u0 

,m 

2 

’pr 

ul+u2 

,m  j 

u2  "  U1 

Pr 

1  u2  +  ul 

Vur u2' 


z  >  0 


(2.33b) 


(2.34a) 


(2.34b) 


z  <  0 

where  sgn[(z-z')]  in  (2.33b)  is  the  numerical  sign  of  (z-z‘). 

With  the  field  transforms  thus  specified,  we  may  employ  (2.10) 
directly  in  order  to  generate  the  desired  integral  representations  for 
the  cartesian  components  of  the  electric  and  magnetic  fields  in  the 
covering  and  substrate  regions,  hence 

k  2r  j-1  .CO 

4®(r,z)  «  -  ;  r(kQar')Hi2)(koar) 


*?(r.z)  -  j;(k0ar')H^>( 

i4n2  ie-iir  v  0 

[2u2n2  -j  -k0(u1z-u2z')  2 

- - - -  e  a  da  > 

L  u0n,  +  u,n~  J 


z  >  0 


(2.35a) 


u2"l  12 

k  rr' 

5!(r,z)  - 


-°-§—  .  0'(knar')H(2)(k  ar) 

j  /i  n 2  J  -ITT  v  cr  v 

i4n^  °°e 


u?ni  -  uinp  k  u  (z+z1)  ~k  up j  z-z*  J  7  2 

-^4 - ^  e  0  2  +  sgn[(z-z')]e  0  2  a2da 

u,nf  +  u,n,  J 

2112  z  >  0  (2.35b) 


iyR  k  f 

*r‘(r,z)  = - 4°  —  i  .  Jv(ko  r')h  (koar) 

i-iir 

ooe 


>  e-ko“2<~'>  *41  ,  2>0 

ul+u2  “2 


.  Ill/  \ 

r(r,z)  = 


iyR  k 


f-2-  J  ,(k  ar')H^(k  ctr) 

2  •  \r  o  '  o 

j  - 1 1. 
ooe 


VU1  kou2<2+2'>  ,  -V212-2'1  (  ad. 
- — t e  +  e  - 


u2+ul 


(2.36a) 


-  ,  z  <  0 


(2.36b) 


where  (2.35a),  (2.35b),  (2.36a),  and  (2.36b)  are  integrated  along  the 
contour  shown  in  Figure  2.1. 

Finally,  with  reference  to  (2.27a)  and  (2.27b),  we  see  that 
(2.35a)  and  (2.35b)  can  be  related  to  u*)(r,z)  by 


«®(r,z)  «  r'  —  p-  n®(r,z) 
r  3r'3z'  r  z 


(2.37) 


2.3.4  ^-Directed  Phased  Arrays  of  Polarization  Current 
As  per  Section  2.3.3,  we  consider  the  cartesian  components 
=  Ez  and  ^  =  Hz  of  the  electric  and  magnetic  fields  in  the  cover¬ 
ing  and  substrate  media.  The  procedure  for  finding  integral  represen¬ 
tations  of  these  cartesian  field  components  is  very  similar  to  that 
followed  above.  Indeed,  the  only  formal  difference  between  the  case 
currently  under  consideration  and  that  of  the  r-directed  polarization 
currents  is  that  the  unit  vector  in  the  representation  of  the  polarization 
current  densities  changes  from  ar  to  a^  ,  with  the  result  that  (A. 3c) 
and  (A. 3d)  [cf.  Appendix  A]  are  replaced,  respectively,  by 
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3z  V  ' 


8  m  /  %  i 

H  V“’z) 


=  r ‘0 '  (k  ar') 

0  V  v  0  ' 


Z=Z  *  + 


Z=2  1  - 


(2.38a) 


$®(a,z‘  +  )  -  4®(a,z'-) 


Wi:on2 


%  Jv(koar,) 


(2.38b) 


Comparing  (2.38a)  with  (A. 3c)  (recalling  v  =  YkQR0)  we  see 


-iyR  k3J  (k  ar 1 )  ->  k3ar'J‘(k  ar‘) 

0  0  V  0  '  0  V  O 


(2.39) 


and  comparison  of  (2.38b)  with  (A. 3d)  shows  that 


k  S  ar‘  R  k* 

-5—9 -  J^(koar') 


R  k^r 

Jv(koar') 


(2.40) 


The  changes  indicated  in  (2.39)  and  (2.40)  manifest  themselves  in  the 

values  of  Bm  and  Be  .  Thus,  (2.32a)  and  (2.32b)  are  replaced  by 

Pi})  Pb 

2 

-  k  ar 

B'  =  — - J’(kar')  ( 

p$  CU2  v  .  0 

yR  k2; 

B*  =  j  (k  ar')  ( 

p*  2n2  v  0 


(2.41a) 


(2.41b) 


For  the  case  of  ^-directed  polarization  currents,  equations  (2.33a), 
(2.33b)  and  (2.34a),  (2.34b)  are  reproduced  except  that  the  subscript  r 
is  replaced  by  the  subscript  $  .  Using  (2.33a)  through  (2.34b)  modified 
for  (^-directed  currents  in  (2.10)  we  have 


f 


rr 
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^(r,z)  = 


yR  C 

0  0  0 


4n 


2  °°e 


-ITT 


Jv(koar‘)H^2)(koar) 


2u0n 


2 
2"2 


2T7 2 

u2nl  +  uln2 


■kn(ul2*u22'> 
e  0  ada  , 


2  >  0 


^e(r>2)  .  Y Voko  ' 


4> 


4n 


2  °°e 


-ITT 


Jv(koar,)Hv2)(koar) 


p  p  ( 

u»nf-u,n9  k  u  (z+2‘)  -ku?\z-z‘\ 

-  '  \--J  e  +  s9nL(2-2‘  )]e  0  2 

u2nl +  uln2 


2  <  0 


.  ni/  \  _  _ u 

V  ’  ~ 

2u 


kZr'  ' 

V  f 


j  _ilT  Jv(koar,)Hv2)(koar) 

cog 


2  ”koU2^2"2 ' ^  a2da 

G 


VU2 


2  >  0 


-k2r ' 


<Ar,z)  =  - -  ^  0‘  (k„ar'  )H^(kAar) 

4>  H  )  .  i 


V  J  -ITT 

oog 
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(2.42c) 
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Upon  comparison  of  (2.42a)  and  (2.42c)  with  (2.27a)  and  (2.27b), 
respectively,  we  see 

*®(r,2)  =  3§.  n*(r,z)  (2.43) 

Finally,  comparison  of  (2.42c)  and  (2.42d)  with  (2.36a)  and  (2.36b), 
respecti vely,  yields 


<(r,2)  =  -v 

<p 


-r' 


iyR  k  3r‘  r 
o  o 


3r  *"(r,z) 
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2.3.5  The  Ancillary  Potential  Integrals  Vg , 


In  terms  of  the  Hertzian  potential  of  Section  2.3.2,  we  can 

p 

write  the  associated  cartesian  component  of  the  electric  field  as 

*z(r,z)  =  (~^"2  kis2^nz  (2.45) 

3  z 

Now,  considering  (2.27a)  and  (2.27b)  we  define  the  scalar,  electric, 
ancillary  potential  Vg  by 


Vp(r,z)  =  (-if-)  ' 
e  4nZ 


-ITT 


Jv(kor‘a)K(2)(kom) 


2u2n2 


2  °°e 

ko^22,“ul2)  a  da 
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(2.46a) 
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(2.46b) 


Having  thus  defined  V  ,  we  see  that 

e 

*®(r,z)  =  r'Mx  +  k^  ol  V  (r,z)  (2.47a) 

z  |3zz  i,Zj  e 

In  a  like  manner,  considering  (2.46a)  and  (2.46b)  along  with  (2.35a) 
and  (2.35b),  we  have 

,2 

Ar, z)  =  r' - V  (r, z)  (2.47b) 

r  3r'3z‘  e 

Consider  now  the  definition  of  the  scalar,  magnetic,  ancillary 


potential  Vm  given  as 
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e  +  e  >  — 
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(2.48b) 


Comparison  of  (2.48a)  and  (2.48b)  with  (2.36a)  and  (2.36b)  gives 


t  (r, z)  =  (iyR  k  )V_ 
r  '  o  o  m 


(2.48c) 


and  a  similar  comparison  of  (2.48a)  and  (2.48b)  with  (2.42c)  and  (2.42d) 
yields 


*;cr.*>  -  -r'  Vm 


(2.48d) 


We  can  conveniently  summarize  the  representation  of  the  total 

p  ft] 

cartesian  field  components  $  and  4  due  to  polari2ation  currents  in 
the  channel  region  in  terms  of  the  potentials  Vg  and  V^,  i.e.,  since 

, e  _  , e  e  .  . e  .  .m  .m  ,  ,rn 

<t>  —  <j>  +  0,  and  4>  =  <S>  +  4>  , 
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CHAPTER  3 

asymptotic  evaluation  of  the  ancillary 

POTENTIAL  V 

e 

3.1  Introduction 

The  developments  in  Chapter  2  have  led  to  a  representation  of  the 
Cartesian  field  components  exterior  to  the  curved  channel  waveguide  in 
terms  of  the  fields  of  the  straight  channel  waveguide  and  tne  ancillary 
electric  and  magnetic  potentials  Vg  and  Vm>  respectively.  In  this 
chapter  we  will  concentrate  upon  determining  an  analytical  expression 
for  Vg  that  is  more  useful  in  practical  applications  than  is  the  formal 
representation  in  equations  (2.46a),  (2.46b).  The  ancillary  potential  \L 

fill 

will  be  discussed  only  to  a  limited  extent  since  the  analytical  procedure 
employed  to  simplify  Vg  can  be  applied  without  modification  to  Vm- 

The  analytical  approach  to  finding  a  more  useful  representation  for 
Vg  and  Vm  begins  with  application  of  the  Debye  expansions  for  Bessel 
and  Hankel  functions  of  large  order  and  argument  to  the  integrands  of 
(2.46a),  (2.46b)  and  (2.48a),  (2.48b).  Once  accomplished,  this  step  is 
seen  to  facilitate  an  asymptotic  analysis  of  the  modified  integrals.  In 
particular,  a  steepest  descent  evaluation  is  pursued  for  Vfi  which 
yields  a  first  order  representation  when  kQRo  »  1. 

In  Chapter  4,  the  asymptotic  expression  derived  for  Vg  is  used  to 
compute  the  attenuation  constant  due  to  continuous  radiation  loss  from 
a  bent  rectangular  waveguide  after  assuming  a  specific  form  for  the  straight 
guide  modes. 


- - jj,-!,  |  -  ' 
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3.2  Modification  of  the  Integral  Representations  of  Vg, 

3.2.1  Changing  the  Contour  of  Integration 

As  we  have  noted  previously,  the  contour  over  which  the  integral 
representations  of  Vg  and  are  integrated  is  that  -illustrated  in 
Figure  2.1.  We  now  consider  the  closed  contour  Fc  shown  in  Figure  3.1 
which  is  composed  of  the  contour  yq  which  lies  along  the  Re  a  axis 

as  well  as  the  contours  y-|>  y2>  Y^  »  Y^  »  YR  •  The  contours  y,  and 

1  2  '3 

Y2  lie  along  the  branch  cuts  associated  with  n-j  and  n^,  respectively, 

and  contours  YD  ,  YD  ,  and  yd  lie  on  the  radius  of  a  semi-circle  of 
K1  k2  k3 

radius  R.  We  will  eventually  consider  T  as  R  -►  whereupcr 
will  coincide  with  our  original  contour  of  integration  for  V  and  » 

For  the  moment,  however  consider  the  integrals 


IG* 


Mv‘  -  uiz) 

F(ot )  e  0  2  1  da  , 


z  >  0 


c 

r  f  k  u  (z+z1)  -k  u  [ z-z 1  ) 

=  i>  |G(a)e  0  +  H(a)e  1  da  , 


z  <  C 


(3.1a) 


(3.1b) 


where  F(a),  G(a),  and  H(a)  are  analytic  everywhere  within  and  on  and 
tend  uniformly  to  zero  as  R-+°°.  Since  there  are  no  singularities  of  the 
integrand  contained  within  r  ,  the  residue  theorem  assures  that  I^=  I^=  0, 
so  (3.1a)  and  (3.1b)  give,  respectively, 

f  (FJ  =  -  i  (F)  -  f  (F)  - 


-  f  (F) 


(3.2a) 


'o  rl  '2 

r  r  «2  K3 

(G,H)  =  -  j  (G,H)  -  j  (G,H)  -  f  (G,H) 
T-,  Y2 


R2  R3 


(3.2b) 
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where  we  hzve  symbolized  the  integrand  of  (3.1a)  by  (F)  and  the  integranc 
of  (3.1b)  by  (G,H). 


If  now  we  let  R  -►  we  may  invoke  Jordan's  lemma  to  justify  tne 


claim  that  both  (F) 

yr,+yr,'V 

l  C  J 


and 


(G,H) 


[ 


vanish  (recall  z‘  <  0). 


Also,  the  uniform  tendency  toward  zero  of  F(a),  G(at),  and  H(a)  together 

with  the  exponential  factors  in  (F)  and  (G,H)  assure  the  convergence  of 

the  remaining  integrals,  whence  (3.2a)  and  (3.2b)  become 
00 


(F)  -  '  (F) 

Y2 


(3.3a) 


00 

r 


(G,H) 


(G,H) 

y2 


(3.3b) 


where  y1  and  are  now  contours  which  "wrap-around"  the  full  extent 
of  their  associated  branch  cuts. 

Notice  that  both  V£  and  Vm  take  the  form  of  (3.3a)  for  z  >  0, 

and  that  both  V  and  V  take  the  form  of  (3.3b)  for  z  <  0.  Thus, 
e  m  — 

we  can  now  transfer  our  attention  from  an  integration  of  the  ancillary 
potentials  Vg  and  along  the  Re  at  axis  to  a  pair  of  integrals 
which,  in  each  case,  follow  the  branch  cuts  of  u-j  and  u ^  in  the  lower 
hal f  of  the  a-plane. 

3.2.2  Application  of  the  Debye  Expansions 

Our  objective  in  this  section  is  to  find  an  asymptotic  representation 
for  the  product  Jv(kQr'a)H^ (kQra)  which  appears  in  the  integrands 
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of  the  integral  representations  of  Vg  and  Vm-  Since  we  are  interested 

in  geometries  for  which  k  R  »  1,  and  since  v  =  vk  R  ,  we  are  certanly 

0  0  oo 

dealing  with  Bessel  and  Hankel  functions  of  large  order.  Likewise,  since 
we  are  interested  in  fields  just  beyond  the  outer  radius  of  the  channel 
guide,  both  r‘  and  r  are  on  the  order  of  Rq;  so,  for  a  of  suffi¬ 
ciently  large  modulus,  the  Bessel  and  Hankel  functions  we  encounter  are 
of  large  argument  as  well.  The  requirement  that  a  be  of  sufficient 
modulus  means  that  there  is  a  region  about  the  origin  of  the  a-plane 

of  radius  within  which  we  cannot  allow  our  integration  contour 

o  o 

to  pass.  This  presents  no  significant  difficulty  since  we  are  free  to 
deform  the  contours  Y-|  and  y2  0nt0  the  improper  Riemann  sheet  (if 
necessary)  for  a  finite  distance  and  hence  to  circumvent  the  forbidden 
region,  whereupon  we  re-enter  the  proper  sheet  and  continue  integrating 
along  the  branch  cuts.  We  will  not  consider  this  subtlety  in  greater 
depth  since  we  will  ultimately  prosecute  a  saddle  point  evaluation  of 
our  integrals  for  which  the  saddle  point  is  well  removed  from  regions 
where  the  large  argument  assumption  is  invalid. 

The  Debye  asymptotic  expansions  for  large  order  and  argument  may 
be  written  in  the  form  (cf.[5],  [9]) 


(2ttv)\(vz)  %  (i  -  z2rVvf(z) 

(2irv)Sv(vz)  %  -2(l-z2)"Vf(z) 
where  0  <  Re  z  <  1 ,  and 

f(z)  =tanh_1(l  -  z2)*5  -  (1-z2) 


(3.4a) 

(3.4b) 

(3.4c) 


and  where  Yy  is  the  Neumann  function  of  order  v  . 
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bince  H^'  =  -  iY  ,  (3.4a)  and  (3.4b)  give 


u  (2~vf 


j2j hit  )~4  j"pvf(z)  _  i  -vf(z) 


(3.5) 


If  we  let  vz  =  k  a(R  +y),  tnen 
o  o 


2  =7°  +  f  > 
0 


(3.6) 


Z„  ‘i  2 

0 


(3.7) 


Recall  that  y  is  the  radial  coordinate  of  the  local  coordinate  system 
situated  at  the  center  of  the  channel  guide  (cf.  Figure  1.2).  In  an 
integrated  optics  environment  we  will  have  n  >  y  >  n^,  thus  inspection 
of  the  contour  reveals  that  Re(^)  <  1  for  all  a  on 

A  Taylor  series  expansion  of  (3.4c)  in  powers  of  y  about  zQ  as 
given  in  (3.7)  yields  the  result 


vf(z)  *  k0R0yf(^)-yk0y(l-^) 

Y 


(3.8) 


Using  (3.8)  in  (3.5)  gives 


H^(vz)  'v  — 

(2ttv) 


n  -  %  (i  *£-> ) 


(3.9) 


to  0(i) 
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With  che  definitions 


x  4  (y2  -  a2) 


(3.10a) 


F(A)  4  yt(~)  =  Y  tanh_1(^-)  -  A 


(3. 1Gb) 


(3.9)  gives  (recall  kQRo  »  1) 


Hv  (koar)  "koAy 

- ^  e  -  ie  sinh  (k  Ay) 

H  2  (k  aR  )  0 

v  o  o 


(3.10c) 


where  we  have  also  used  r  =  R  +  y. 

o 

Using  (3.4a)  with  (3.6),  (3.7),  and  (3.8),  we  find 


Jv(koarH 


2  2  -h 

n  +/-)  j 


(2nv) 


k0*y(l-^)  -Yk0V(Y> 

!  Y  e  0  0  Y  (3.11) 


to  O(^).  Employing  (3.11)  we  find  directly 


J  ,(k  ar) 
v  o 


v'  o  o' 


Also,  using  (3.9)  and  (3.11)  we  have 


/?1  l  fi  -2Yk„Rrtf(a) 

Jv<ko“VHv  <koaRo>  '  STrT j?  e  Y  +  1 

0  0/ 


(3.12) 


(3.13) 


We  now  have  each  of  the  elemental  results  necessary  for  finding 
an  asymptotic  representation  for  Jv(k0ar*  )H^ (kQar) .  From  (3.10c)  and 
(3.12)  we  have 
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y'k  x 


J  (k  ar')H^(k  ar)  a,  J  (k  aR  )H^(k  aR  )e  0 


v  o  v  o 


V  o  o  v  'oo 


-k  Xy  -2k  R  F(X) 

?  0  -  i  e  00  sinh(kQXy) 


Finally,  substitution  of  (3.13)  in  (3.14)  yields  our  objective, 

i2)  i  fi  -kX(y-y') 


J  (k  ar  )H^  ^(k  ar)  'v  -p- -  j  e 
v  o  '  v  v  0  TTk  R„  )  X 

o  o  / 


1  V(y+y,)  ~2koRoF(X) 


+  2X  e 


-4koRoF{A) 

where  terms  of  0  (e  )  have  been  neglected. 


(3.14) 


(3.15) 


3.2.3  The  Modified  Integral  Representations  of  Vg, 

Using  results  (3.3a),  (3.3b)  and  (3.15)  in  (2.46a)  and  (2.46b)  as 

well  as  (2.48a)  and  (2.48b),  the  ancillary  potentials  Vg  and  seen 

to  be  separable  into  primary  and  secondary  components ;  i.e.,  we  may 

write  V  =  Vp  +  vf  ,  V  =  +  vf.  where 

e  e  e  m  m  m 
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2  '  4-rrk  R 
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2  Y 


-k0X(y-y‘ ) 
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(3.16a) 


;-}(3. 16b) 


(3.16c) 


Vm(y  ,y  *  ,z,z' )  =  (jdp) 
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( 3. 16d) 
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For  convenience  we  have  omitted  detailed  representation  of  those 
portions  of  the  integrands  which  remain  unchanged  as  the  result  of  apply¬ 
ing  the  Debye  expansions.  Classification  of  terms  as  primary  or  secondary 

components  is  made  on  the  basis  of  the  fact  that  in  the  limit  as  R  co, 

0 

RqV^  and  RqV P  approach  the  forms  assumed  by  the  ancillary  potentials 

in  the  straight  guide  case.  Because  of  the  presence  of  the  factor 
-2k  R  F(A) 

e  ,  choice  of  the  proper  branch  of  A  will  in  general  yieid 

VS«VP  and  Vs  <r<  \/P;  indeed,  both  R  Vs  and  R  Vs  vanish  as  R  ->■  «  . 
ee  mm  oe  om  o 

Despite  the  fact  that  the  secondary  ancillary  potentials  are  much 

s  s 

smaller  in  magnitude  than  their  primary  counterparts ,  Vg  and  V„  are 

the  terms  of  interest  for  our  purposes  here.  This  follows  from  the 

fact  that,  in  addition  to  their  relatively  small  magnitudes  when  compared 

to  the  primary  potentials,  both  and  exhibit  an  exponentially 

increasing  character  for  increasing  y,  and  each  stands  in  quadrature 

to  its  associated  primary  potential.  These  features  of  and 

allow  an  association  between  these  secondary  potentials  and  the  electric 

and  magnetic  fields  exterior  to  the  bent  channelguide  which  have  been 

reflected  from  the  caustic;  i.e.,  the  electric  and  magnetic  fields 

s  s 

derivable  from  the  secondary,  ancillary  potentials  Vg  and  Vm  are 
the  fields  dE  and  6H,  respectively. 

We  will  henceforth  focus  attention  on  the  ancillary  potentials  Vg 
of  (3.16b)  and  of  (3.16d)  in  order  to  characterize  the  radiation 
from  the  bent  channel  guide.  Because  of  the  linearity  of  all  operators 
involved,  equation  (2.49)  still  applies  and  the  components  produced  are 
the  z-components  of  <SE  and  bH. 
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3.3  The  Asymptotic  Form  of  \IQ 
3.3.1  A  Change  of  Variables 

To  assure  convergence  of  our  fundamental  integral  representations  of 

2  p  h 

Vg  and  Vm ,  we  have  chosen  to  employ  the  branches  of  u-  =  (a  -  n-) 

2  ?  ^ 

and  u0  =  (a  -  n£)  for  which  Re  u^  >0,  Re  u^  >  0.  It  is  in  connec¬ 
tion  with  these  branches  of  u-j  and  u^  that  the  branch  cuts  at  n.  and 
n^  arise  in  the  a-pl ane.  The  contours  and  y^  follow  the  o ranch 
cuts  at  n-j  and  n^,  respectively,  and  it  follows  from  the  same  analysis 
that  specifies  the  branch  cuts  (not  discussed  herein)  that  u.  is  purely 
imaginary  along  y^  and  u^  is  purely  imaginary  along  X  .  Indeed,  if 
we  consider  y^  to  be  composed  of  segments  y^,  g^ >  Y2u  (i-e., 

y2  =  y2 i  +  92  +  y2u^  is  Possible  t0  show  that  ar9  u2  =  "  |  on  Y->jo 

and  arg  u^  =  ^  on  y^u  (cf.  Figure  3.2).  The  segment  g^  is  a  circular 

segment  centered  on  n^  whose  radius  we  allow  to  vanish,  whereupon  we  find 

no  contribution  from  g^  to  the  integrals  over  y^  in  (3.16b)  and  (3.1 6d) . 

Analogous  statements  follow  for  the  argument  of  u-j  on  Y^ ,  with 

y-j  =  y-j^  +  g-j  +  y-ju,  and  for  the  contribution  from  g^  to  the  integrals 

over  y.  in  (3.16b)  and  (3.16d). 

2  2  *5 

Analysis  shows  that  we  may  select  a  branch  of  the  function  {r\%  -  a  ) 

2  2  *5 

consistent  with  the  extant  branch  cut  at  n2  such  that  Im(n2  -  a  )  <  0, 
thus  we  may  write 


u2  =  1  s2 


(3.17a) 


where 


oo  2  2  ^ 

s2  =  (n2  "  a  '  ’  Im(n2  "  a  )  <  0 


(3.17b) 
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Since  is  purely  imaginary  along  and  Y2u»  then  S2  must  be 

purely  real  along  these  same  segments.  Indeed,  the  nature  of  u 2  along 
and  Y2U  along  with  (3.17a)  requires  Sg  real  and  s?  <  0  on 
Y2£  and  S2  real  and  S2  >  0  along  Y2u  (of-  figure  3.2).  Thus,  a  change 
of  variables  from  a  to  S2  in  (3.16b)  and  (3.16d)  using  (3.17a)  will 
convert  integrals  along  Y2  in  the  a-plane  to  integrals  along  the  full 
extent  of  the  real  axis  (-»  to  +°°)  in  the  s^-plane.  In  the  same  way, 
one  may  write 


U1  =  1  S1 


where 


si  -  (ni 


Im(n, 


2  ** 

a)  <  0 


(3.18a) 


(3.18b) 


as  a  change  of  variables  from  a  to  s^  in  (3.16b)  and  (3.16d)  to  convert 
integrals  over  y-j  in  the  a-plane  to  integrals  over  the  full  extent  of 
the  real  axis  in  the  s-j-plane. 

For  the  purpose  of  illustration,  we  will  consider  that  portion  of 
V*  obtained  from  integration  over  the  contour  Y2>  while  noting  that  an 
analogous  procedure  may  be  performed  for  that  portiom  of  v^  obtained 
from  integration  over  y-j*  At  the  end  of  this  chapter  we  will  see,  in 
fact,  that  the  contribution  to  V*  due  to  the  integration  over  yt  can 
be  neglected  because  of  its  insignificant  amplitude  relative  to  the  con¬ 
tribution  arising  from  the  integration  over  Yg 

=  -ia(at  -  n£) 

ada 


From  (3.17a)  and  (3.17b)  we  have  ^  -  J-'-Z  2 


-h 


so 


ds2  =  -i 


(3.19) 
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3.3.2  Evaluation  of  V' 


Since  we  assume  kQRo  »  1,  the  form  of  Vg  given  in  (3.22a)  lends  itself 
readily  to  a  steepest  descent  type  of  asymptotic  evaluation  (cf.  [7]). 


Proceeding  as  outlined  in  Appends  B  we  find 

/  _2  l3'2 


Vg(y.y'.z.z')  •x. 


16Tr2nf k  R 
2  o  o 


A„k  R 
2  o  o 


Voq2(0) 


,(y+y' ) 


2  /  2  2  7% 

^2 ^ i / 


koz  - 


n2 /  2  2,H 

^2^2  ”  ^1 


(3.23) 


A  similar  result  is  obtained  for  the  component  of  derived  by 
integration  over  .  However,  consider  the  series  representation  of 
q-j ( 0)  and  q2(0)  derived  by  expanding  their  logarithmic  terms  [cf.  (3.22d)J 


qi(°)  ~  ~2y  ^(2n+l)  ^  y  ^ 


A .  2n+l 


(3.24) 


where  X.  n2  ^ 

(~)  =  O  -  4“)  ’  i  =  1,2  • 

Y  Y 

n  2  A-j 

For  q-j  (0) ,  (—-)  «  1,  so  (-—■)  ~  1  , 

n?  2  A„ 

whereas  for  q2(0),  ~  1,  so  0  <  (■—)  «  1.  For  typical  integrated 

optics  material  parameters  then,  (3.24)  gives  [q-j  (0)  |  »  |q2(0) ,  so  that 
knRoql(0)  kAq2(0) 

e  «  e  c  and  we  are  justified  in  neglecting  the  contri¬ 


bution  from  the  integration  over  contour  y^. 
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CHAPTER  4 

RADIATION  LOSS  FROM  A  CURVED  RECTANGULAR 
CHANNEL  GUIDE  FAR  FROM  CUTOFF 

4. 1  Introduction 

As  discussed  briefly  in  Chapter  1,  we  will  know  the  continuous 
radiation  loss  from  a  section  of  curved  dielectric  channel  waveguide 
once  we  have  found  the  normalized  attenuation  constant  a  as  given  by  (1.2). 
Evaluation  of  a  requires  computation  of  the  quantities  P  and  c  in 
(1.3)  and  (1.4),  respectively.  Our  efforts  in  Chapters  2  and  3  have 
been  to  derive  an  expression  from  which  we  can  determine  the  caustic- 
reflected  fields  61  and  6H  in  terms  of  the  straight  channel  fields. 
Equations  (2.49)  and  (3.23)  allow  us  now  to  do  so.  Thus,  we  are  left 
with  the  task  of  specifying  the  straight  channel  fields  EQ,  Ho  for  any 
case  of  practical  interest  and  using  these  with  (2.49),  (3.23),  (1.2),  and 
(1.3)  in  order  to  find  a. 

The  general  problem  of  determining  the  fields  within  a  straight, 
rectargular  dielectric  channel  waveguide  embedded  in  a  substrate  has  not 
been  r.gorously  solved  analytically.  However,  an  approximate  analytical 
method  has  been  developed  by  Marcatili  (cf.  [2])  which  yields  straight 
channel  field  expressions  which  are  valid  for  well  confined  channel  modes, 
i.e.,  channel  modes  that  are  far  from  cutoff.  We  will  employ  Marcatili’s 
results  and  thus  generate  an  expression  for  the  radiation  from  a  curved 
dielectric  channel  waveguide  with  associated  straight  channel  modes  that 
are  far  from  cutoff.  In  particular,  we  consider  the  straight,  rectangular 
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dielectric  channel  waveguide  configuration  of  Fig.  4.1  which  illustrates 
the  general  disposition  of  material  regions  to  which  Marcatili's 
analytical  procedure  can  be  applied.  The  material  parameters  assigned 
to  tne  various  regions  are  consistent  with  those  or  the  assoc ucec 
curved  channel  guide  problem  illustrated  ir,  Figure  1.2.  The  assumption 
of  propagating  modes  that  are  far  from  cutoff  and  therefore  highly  con¬ 
fined  within  the  channel  region  is  extended  by  the  stipulation  that  the 
fields  in  the  shaded  regions  of  Figure  4.1  are  negligible,  hence  no 
material  parameters  need  be  specified  there. 

It  should  also  be  noted  that  the  origin  of  the  (x,w,s)  coordinate 
system  is  off-set  from  the  local  cartesian  system  (x.y.s)  ir,  Figure  1.1. 
To  convert  to  the  (x,w,s)  system  from  the  (x,y,s)  system  we  need  only 


employ  the  transformation 


y  =  w  -  g 


(4.1) 


The  structure  of  Figure  4.1  can  support  a  set  of  modes  that  are 


predominantly  x-polarized  and  are  designated  modes.  In  addition, 
such  a  structure  can  also  support  a  set  of  modes  that  are  predominantly 
y-polarized  and  are  designated  E^  modes.  Although  such  a  structure 
can  support  both  modes,  we  consider  only  radiation  due  to  the  mooes 
whose  field  distributions  in  the  channel  guide  are  given  by  Marcuse 


(cf.  [11]) 


Eqs  =  A  cos  ^(x  +  £)cos  kw(w  +  n) 


2  2  2 
c  =  iA  n  ko  -  kx 
ox  S  kx 


sin  kx(x  +  £)cos  kw(w  +  n) 


(4.2a) 


(4.2b) 


E  -  0 
ow 


(4.2c) 


" 1 . >!■ 
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Hos  =  _A(r— )  (j^O (r~)  sin  kx(x+i)sir.  k^w  +  n) 
'  0  x  K 


(4. 2d) 


H  -  0 
ox 


n  n 
o 


ow 


iA(j—  )  sin  k  (x  +  Ocos  k  (w  +  r,) 


x  o 


(4.2e) 
( 4 . 2  f ) 


where  kQ  =  w/mqcq  is  the  free-space  phase  constant,  n  is  the  phase 
constant  of  the  fields  in  the  channel  (propagation  of  the  form  e"‘'wt~"' 

is  the  free-space  wave  impedance,  "he  para- 


is  assumed),  and  c  = 

U 

meters  kx  and  k  f  are  separation  constants  related  by 


2  2  .2  ,2  ,  i  2 

n  xo  -  ^  =  kx  +  ky 


(4.3) 


and  are  solutions  of  the  eigenvalue  equations 

2  2 

n  k  (n-jYo  fioYo) 
tan  k  d  -  -A-y  - — -  -  -  -  - 
(n1n2kx  -  n  y2Y3/ 

2k  y, 

tan  k  b  =  — ~ — — 
w  -l 


~  Y/i ) 

w  1 4 


\4.4a) 

(4.4b) 


The  parameters  Yi  ace  given  by 
0 


y2  *  t(n2-  n|)k \  -  kllH 
n  *  t(a2  -  o2X  - 


(4.5 a) 
(4.5b) 
(4.5c) 


(3.18a) 
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The  parameters  y.  are  transverse  decay  constants  of  the  fields  in 

the  material  regions  below,  above,  ana  on  either  side  of  the  region  of 

refractive  index  n  (the  channel  guide),  respectively;  e,g.,  the  fields  in 

Y?(x+d) 

the  material  region  below  the  channel  guide  region  decay  like  e  , 


x<  -d.  It  follows  that  each  of  the  parameters  y.  satisfies  y-  >  0. 
~  Jo 

Once  kx  and  kw  are  known,  the  phase  parameters  £,  and  r,  are 


determined  from 


tan  k  £  =  -  (—)  — 
x  n  '  y  j 


tan  k  n  = 
w 


(4.6a) 

(4.6b) 


Also,  notice  that  once  k^  and  k w  are  known,  (4.2)  may  be  solved  to 
yield  3- 


4.2  Expressions  for  P  and  C 


Application  of  the  assumption  of  modes  far  from  cutoff  and  employing 


the  local  coordinate  representations  for  the  field  components  allows 
(1.3)  and  (1.4)  to  be  rewritten,  respectively,  as 
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o  b 
P  =  2  .  dx  f 


■  dx  dy(£  H  ) 
J  j  ox  ow' 


■\U 


-d  o 


x  6H  -  6E  x  H"  1  •  a  dx 


o  f  w 


(4.7) 


(4.8) 


As  discussed  in  detail  in  Appendix  C,  the  straight  channel  field 
expressions  (4.2)  nay  be  used  with  the  ancillary  potential  from  (3.23) 
in  (4.7).  and  (4.8)  to  evaluate  P  and  C  as 


P  =  -2A 


2  2  ?  2 

0  rnr  -  kS  k  rr  ,  K 

2  _ S _ x  j  _ Q_  +i  +  i  \ 

L  k2  J  eco  '  2  LXH2  V 


(4.9) 


c  ”  'a2(:1C2  s9nLcos(kwb)jsgnCk^  -  Y^j 


■WAo  -  -  "!>1 


rx.  x 
w  o 

16m 


2  2  2 

1 677^3  Vx2?0 


[V^O  -  *x>  +  "V<] 


k  R  A, 
o  o  2 


(k  A2b)/2  koRoq2(0) 
e  e 


(4.10) 


where  the  constants  Lx,  Lw>  C^,  and  C2  are  given  in  Appendix  C  as  equations 
(C.5a),  (C.5b) ,  (C.22) ,  and  (C.25),  respectively. 


4.3  The  Evaluation  of  a 

The  normalized  attenuation  constant  is  given  in  (1.2).  Before  evaluating 
a  in  terms  of  (4.9)  and  (4.10),  it  is  advantageous  to  normalize  Cy  C2>  anc 
A  -  (4  +  L  )(-«-  +  L  )  to  k  .  Doing  this  we  have 

G  t  X  t  W  0 

C>  ■  .“27 zrj?  {A2C59n[coS(k0l>)vw]s9n[v2  - 
K0r2(n  -n2)  { 
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where 
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=  1~4  !-)  fx,(sgn[cos(k  b)v  )]sgn[v®  - 

N  ["T‘j  1 

'kob)'2/2  _  '(kob'X2/2)  *  r  (s9n[cos(k0b)vw3 


sgn[v 


2  .  j.2-,  e(kob,X2/2  +  e‘'<ob,X2/2)‘j  [sgn[cos(k0<Dv, 


(n2r2(k0“)  ’  "i1  . 

sgn[n2n|v2  -  "  ~2r33  \  jJTJ3F~ 


1 _ , 

,42  4  2\i 

(n  r2  +  n,vx) 


(4.16b) 


(4.16c) 


(4.17b) 


In  order  to  highlight  the  dependence  of  this  cumbersome  result  on 
the  radius  of  curvature  of  the  channel  we  simplify  (4.16a)  through  (4.16c) 


by  writing 


(k  R  ) 
'  o  o 


Ci__  AW0* 

o  '4/2 


(4.18) 


where  (the  subscript  i  designating  that  this  result  is  associated 
with  an  "inhomogeneous"  medium  about  the  channel  guide)  is  a  para.„ater 
dependent  upon  the  refractive  indices  of  the  channel  guide  and  its  surround- 
ing  media,  the  dimensions  of  the  channel  guide,  and  the  specific  Ep^  mode 
propagating  in  the  channel. 

Figure  4.2  is  a  plot  of  the  normalized  attenuation  coefficient 

versus  k  R  for  k  d  =  5  and  kb  =  10.  The  refractive  index  of 
oo  o  o 

the  channel  guide  is  assigned  to  be  4.04,  that  of  the  substrate  4.00, 
and  that  of  the  covering  medium  1.00. 

4.4  Conclusions 

As  mentioned  in  Secion  1.1,  there  are  no  analytical  reports  in  the 
literature  which  address  a  problem  similar  enough  to  that  considered 
here  to  allow  comparisons  with  the  force  to  effectively  confirm  or 
reject  the  principal  result  derived  in  this  report,  i.e.,  the  form  of 
the  normalized  attenuation  coefficient  a  presented  in  Section  4.4. 

To  this  extent,  then,  the  expression  for  a  in  (4.16)  is  a  new  result. 

To  appreciate  that  the  form  of  a  constitutes  a  unique  as  well  as  a 
new  result,  we  consider  the  form  of  the  normalized  attenuation  coefficient 
a1  (primes  do  not  indicate  differentiation  here)  for  a  bent,  rectangular 
dielectric  waveguide  immersed  in  a  homogeneous  medium  of  lower  refractive 
index  than  that  of  the  guide  (cf.  [5],  Section  9.8),  i.e. 
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a 


W 


W(0) 


(4.19) 


where  and  Q(0),  like  and  q  (0)  in  (4.41),  are  parameters 
depending  upon  material  characteristics  and  the  specific  mode  assumed  to 
be  propagating  in  the  guide. 

Upon  comparison  of  (4.19)  with  (4.18),  two  conclusions  are  immediate. 
First,  even  though  Q(0)  and  q2(0)  are  not  necessarily  equal,  the  form 
of  the  exponential  dependence  of  a  and  a'  on  kQRQ  is  the  same. 

Second,  the  algebraic  dependence  of  a  and  a'  on  R  is  not  the  same, 

1  3/2  ,3s  0 

being  for  a  and  ( k~ -g-)  for  a'.  Certainly,  then,  fer  a 

o  o  o  o 

pair  of  bent,  rectangular  dielectric  waveguides,  one  in  a  dense  substrate 
below  a  material  interface  and  the  other  in  a  homogeneous  medium,  when 
q2(0)  =  Q(0),  one  expects  the  rate  of  decrease  of  attenuation  with 
increasing  kQRo  to  be  greater  in  the  case  of  the  guide  in  the  substrate 
below  a  material  interface.  In  order  to  understand  how  the  overall  atten¬ 
uation  behaves  in  such  situations  we  consider  Figure  4.3,  which  is  a  plot 
(versus  kQRQ)  of  the  values  of  the  normalized  attenuation  constant  for 
each  of  three  bent  waveguide  configurations. 

In  case  I,  the  guide  has  a  refractive  index  n  =  4.04,  dimensions  of 

k  d  =  3,  kb=10,  and  is  embedded  in  a  homogeneous  medium  of  refractive 

oo 


index  n^  =  4.00.  In  case  II,  the  guide  has  a  refractive  index  n  =  4.028, 

dimensions  of  kQd  =  5,  kQb  =  10,  and  is  likewise  embedded  in  a  homogeneous 

medium  of  refractive  index  n-j  =  4.00.  Case  III  has  a  guide  of  refractive 

index  n  =  4.04  and  dimensions  k  d  =5,  k  b  =  10  embedded  in  a  sub- 

o  o 

strate  of  refractive  index  n2  =  4.00  below  a  material  region  refractive 


-  - 
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index  n-j  =  1.00.  The  dimensions  in  case  I  were  altered  relative  to  tr.ose 
of  case  III  in  order  to  provide  q2(0)  =  Q j ( 0) .  Likewise,  the  refractive 
index  of  the  guide  in  case  II  was  altered  relative  to  that  of  the  guide 
in  case  III  to  again  assure  q^(0)  =  QTT(Q).  Thus,  the  exponential  varia¬ 
tion  of  the  normalized  attenuation  coefficient  is  the  same  for  the  guides 
in  each  of  cases  I,  II,  and  III. 

From  Figure  4.3  we  see  that  in  addition  to  a  more  rapid  decrease  in 
attenuation  with  increasing  koRQ,  The  guide  in  a  dense  substrate  below 
a  material  interface  actually  exhibits  a  lesser  absolute  attenuation  than 
do  the  guides  surrounded  by  homogeneous  media.  Indeed,  while  the  differ¬ 
ence  in  values  of  normalized  attenuation  constant  between  cases  I  or  II 
and  III  is  less  than  an  order  of  magnitude  for  kQRo  =  500,  it  lies 

between  one  and  two  orders  of  magnitude  for  k  R  =  5000. 

oo 

The  divergent  tendency  of  the  curve  for  case  III  relative  to  those 

of  cases  I  and  II  results  from  the  distinctly  different  algebraic 

variation  with  k  R  in  case  III  relative  to  cases  I  and  II.  However, 
o  o 

this  divergent  tendency  is  ultimately  overcome  by  the  identical  exponen¬ 
tial  decay  terms  present  in  each  case.  This  fact  is  illustrated  in 
Figure  4.4,  which  is  a  plot  of  the  normalized  attenuation  constants  for 
cases  I  and  III  over  extended  ranges  in  both  the  abscissa  and  the  ordinate. 

Recalling  the  inverse  three-halves  power  algebraic  dependence  of  a 
on  kQRQ ,  we  may  interpret  this  behavior  as  a  kind  of  lateral-wave 
phenomenon  due  to  the  presence  of  the  material  interface.  This  interpre¬ 
tation  implies  that  the  direct  radiation  which  accounts  for  losses  of  a 
curved  guide  in  a  homogeneous  medium  is,  in  the  case  of  the  guide  in  a 
dense  substrate  below  a  material  interface,  effectively  cancelled  by  an 
image  contribution  from  the  material  interface. 
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Appendix  A 


Determination  and  Application 

of  Source  ana  Boundary 
"6  ""01 

Conditions  for  $  , 


A . 1 .  Introduction 

Inspection  of  (2.30a)  through  (2.30c)  reveals  that  we  nave  eight 
unknown  transform  coefficients,  hence  we  expect  to  be  able  to  specify 
eight  independent  equations  relating  these  coefficients  as  determined 
by  boundary  and  source  conditions  on  tne  field  components.  We  begin  by 
considering  equations  (1.9a)  through  (1.9b)  which  relate  the  transforms 
of  tne  fields  transverse  to  the  z-direction  to  the  transforms  of  the 
cartesian  field  components.  Direct  integration  of  (1.9a)  through  (1.9b) 
after  dividing  both  sides  by  a  and  multiplying  through  by  0  (k^ar)  gives 
the  following  results: 

«x> 

H.(r,z)  =  icQ  j  ?lJ(a,z)J^(koar)da 
0 


-  i  Y  R - 


da 


-tnYR0 
EJr.z)  =  -Vs- 


ou 


0 

5®(o,z)J^(k0or)da 


(A. la) 


(A. lb) 


(A. 1 c) 


'int,2  ^(o^z)J^(k0^r)cta 


’0  0 


V  . 

Hr(r,2)  =  j  *J(«,2)Jv(k0ar)% 

0 


+  r  ’  it  *r(a’zK( Vr)da 

o  J 


(A. Id) 


A. 2  Boundary  Conditions  at  the  Dielectric  Interface 

Since  E  ,  and  E„  must  be  continuous  at  z  =  0,  inspection  of 
r 

(A. la)  and  (A. lb)  yields 


-m 


^!(a,0+)  =  ^(a,0-) 


(A. 2a) 


_3 

3Z 


2=0+ 


3 

3Z 


Z=0- 


(A.2b) 


Likewise,  continuity  of  H?  and  Hf  at  z  =  0  and  inspection  of  (A.lc)  and 
(A. Id)  give 


n?  «^(a,0+)  -  n^  $^(a, 0-) 


(A. 2c) 


A . 3  Boundary  and  Source  Conditions  at  the  Source  Plane 

Tne  continuity  of  and  at  z  =  z‘  and  inspection  of  (A. la)  and 
(A. lb)  yields 


«“(a,z’+)  =  ^(a  ,Z‘~) 


(A. 3a) 


if '<<«•*>' 


■51  »'<«•"). 


(A. 3b) 


‘z  =  z!  + 


1  z  =  z'  ■ 


Using  the  fact  that  the  magnetic  intensity  H  must  be  discontinuous 
as  the  result  of  surface  currents,  while  the  electric  flux  density  D  must 
be  discontinuous  at  points  where  surface  charge  accumulations  exist,  we 
find  the  additional  expressions 


3  ~m,  , 

3Z  vrAa’  1 


si 


=  -ivkoJv(Vr')  (A. 3c) 


>z=z'+ 


‘ z=z*  - 


k^r  ar 1 

*J(a,z'+}  -  «*(a,Z*-)  =  -9-f -  J'v(kor'a) 

in^ 


(A. 3d) 


Equations  (A. 2a)  through  (A. 2d)  and  equations  (A. 3a)  through  (A. 3b) 
constitute  the  independent  set  of  eight  equations  we  need  in  order  to 
solve  for  the  transform  coefficients.  Indeed,  using  (2.30a)  and  (2.30b) 
with  (A.2a)  through  (A. 2b)  we  find 


„  2u  0 

jn  —  2  Bm 

r  u^+u2  p r 


(A. 4a) 


nr  “  u2+u]  Pr 


(A. 4b) 


2^n\ 

r  u2nl+uln2  P 


(A.4c) 


=  u2nru1n2  Re 

u2nl+uln2  P 


(A.4d) 


Likewise,  using  (2.30b)  and  (2.30c)  witn  (A. 3a)  through  (a. 3d),  we  he 


m  f  U0-U,  K  u0z 1  -k  u^z'-)  „ 

Dm  =  e  0  2  +  o  2  Bm 

r  j  u^u2  |  pr 


2  2 
.  fu0n, -u,n, 
ne  =  -LI _ lJL 

ur  i  2  2  e 

(  u2n-j+u^n2 


k  u0z'  -k  u0z‘ 


pr  2u. 


k^c  ar1 


J  (Ur1) 

v  O 


J* (krtr'a) 

v  0 


(A. 3a) 


(A. 5b) 


(A. 5c) 


(A.5d) 


Since  (A. 4a)  through  (A.4d)  and  (A. 5a)  and  (A. 5b)  are  given  in  terms  of 
either  B™r  or  B^,  and  since  these  are  given  explicitly  by  (A. 5c)  and  (A.5d), 

m  p 

respectively,  the  transform  coefficients  $r  and  <?r  are  now  known. 

0  m 

A. 4  Integral  Representations  for  $  ,  » 

Upon  substitution  of  (A. 4a)  through  (A.4d)  and  (A. 5a)  and  (A. 5b) 
into  (2.30a)  through  (2.30c)  we  find 


2u?n2  -k  (u^z-upz1) 

•!(a,2)  =  B«  -|-2— y  e  0  1  2  ,  Z  >  0 

M  u2n^+u1n2 


(A. 6a) 


-m,  _„m  2u2  -“o'VV)  ,  >n 

4r(“>2)  '  Bpr  ^  e  •  2  i  0 


(A. 6b) 


7e,  ne  fu2nrul"2  kou2(2*2'> 

r  pr  (u2n‘+Uln‘ 

-k  u2;2-z*|) 

+  sgnL(z-z‘ )je  ?»  z  <  0 


(A. 6c) 


-r»,  .n  fu2'“i  V2(W!'>  „  •kou2|2-2',J  ,  ... 

“pr  '  +e  t,2jO(A.6o) 


where  sgn l ( -z-z 1 ) 3  in  (A. 6c)  is  the  numberical  sign  of  (z-z1). 

With  the  field  transforms  thus  specified,  we  may  employ  (2.10) 
directly  in  order  to  generate  the  required  integral  representations. 


5r‘r-2>  =  Jv<V')Hl2,(ko“r> 

i4n2  °°e 


2u2n2  -ko(u1z-u2zl)  2 
- -  Q  a  act ,  ,Z  >  U 


u2nl+uln2 


(A. 7a) 


m  i  f 

*>.z).  Jv(ko“r'>VVr> 

coQ 


J^2_  /W2-2')  2>0 

u^Ug  u2 


(A. 7b) 


Appendix  B 


The  Steepest  Descent  Evaluation  of  V' 


To  find  the  saddle-points  of  q2(s2)  1n  s^-plane  we  must  evaluate 
dq2 

,  and  upon  doing  this  we  find 


-[^V] 

Ln,-s0  J 


(B.l) 


From  (B.l  )  we  conclude  that  s^  =  0  is  a  saddle-point  for  q^Sg).  Further 
analysis  shows  that  for  Im  n2  «  1  (as  we  assume),  the  path  of  steepest 
descent  is  essentially  along  the  Re  s2  axis. 

Differentiation  of  (B.l)  provides 


(B.  2) 


Notice  that  since  A  is  an  even  function  of  s2,  then  q(s2)  is  an  even 

3 

function  of  s2,  so  the  coefficient  of  the  s2  term  in  a  Taylor  series 
must  vanish,  thus 


(B.3) 


A  Taylor  series  representation  of  q(s2)  through  the  first  four  terms 
about  s2  =  0  is  thus,  using  (B.l),  (B.2),  and  (B.3), 


(B.4a) 


B-2 


^2^2^  = 


q2(o) 


A2 


where 


q2(0)  =  2  |a2  +  ^  log 


y-A, 

_ C 

Y+Xr 


(B .4b) 


Consider  now  y  for  z<0  as  given  by  (3.22c).  If  we  let  u9  =  is.,. 


then 


Ye(y,y\z,z’is2)  = 


ek0x(y+y  )  pS2n2_u  „2  ik0(z+z’)s2 

A  {  ~  2  2  e 

(_  is2n^+u-jn2 


+  e 


-ikQ  iz-z'i  s2 


(B.5) 


Consideration  of  (B.5)  for  s2  =  0  yields 


y  (y,y‘ ,z,z' ;o)  =  0 


(B.6) 


so  we  must  evaluate 


32Y, 


V 

s2 


s2=0 


in  order  to  generate  the  leading  term  in  the  asymptotic  expansion  of 
V^.  Doing  this,  we  find 


— 


A.) 

4 


-  4?- 


k0(y+y*)tTr2-n|> 


[V  - 


2,  2  2,1/2 

n2Vn2“nl 1 


.z  - 


2/2  2,1/2 

n2'n2“nl 


(B.7) 


The  leading  term  of  the  asymptotic  expansion  of  Vg  has  the  form 


(cf.  [10]) 


-5  k  R  ql^  r32Ye 

Vg( y,y',z,z')' - f -  e  *  [  ~ T 

6  8™2koRo  L  S2 


-kRn/X2;s^  o 
e  0  2  s2ds. 


(B.8) 


Using 


2  2  3/2 

.-■•My  . 


and  (8.7)  in  (B.8), 


Vg(y.y‘ »z,z‘) 


^2"2koRo  V2*0"0 


k0A2(y+y' ) 


_  [y  -  VWwl  [ 

1°  ni,"!-"?) v  J  L 


V  '  n2,  2_  2vl/2j  (6.9) 

n2kn2”nl 


which  is  (3.23)  ih  Section  3.3.2. 


0-1 


The  Evaluation  of  P  and  C 


C.l .  The  Evaluation  of  P 

Application  of  (1.3)  to  the  channel  guide  cross  section  in 


Figure  4. 1  gives 


o  b 


P  =  2  ,  dx  i  dy(Ert  H  ) 

J  ,  J  J  ox  ow' 
-d  o 


(C.l) 


Our  consideration  of  modes  far  from  cutoff  manifests  itself  in  (C.l)  in 
that  the  integration  can  be  taken  over  the  guide  cross  section  only. 
Substitution  of  (4.2b)  and  (4.2f)  in  (c.l)  yields 

9  9  9  9  0  b 

9  n KrK  M  [9  f  9 

P  =-2A^  — sin^Ck  (x+0]dx  J  cos^[k  (w+r.) ]dw  (C.2) 
k  -d  0 

x  0 

f  2 

Expansion  of  the  integrand  in  j  sin  [k  (xH)jdx  and  use  of  (D.7a),  (D.7b), 

-d  x 

(D.9a),  and  (D.9b)  [cf.  Appendix  D]  gives 


0  p 

(9  fi  n2  f  "lY, 

j  .  S1n  kx(x+5)dx  =  2  +  T  /  4  2  4 

-d  (n  Y3+n1 


(C.3) 


Similarly,  expansion  of  the  integrand  of  cos  [k  (w+n)]dw  and  use  of 

w 

(D.4a),  (D.4b) ,  (D.6a),  and  (D.6b)  gives 


cos2kK(w*iudw  ■  \  +  k2(n2.n2) 


(C.4) 


0 


C-2 


Making  the  definitions 


f  nlY3  n2Y2  1 

r  4  2.  4.2  ,  4  2.  4,2,  ( 

[( n  Y3+n1kx  (n  Y2+n21V  > 


(C.5a) 


y4 


,2/  2  2, 

kQ(n  -n2) 


and  substituting  (C.3)  and  (C.4)  into  (C.2)  gives 


(C. 5b) 


(C.6) 


C-2.  The  Evaluation  of  C 

C.2.1  An  Explicit  Integral  Representation  of  C 

22  2  1/2 

For  the  circumstance  of  modes  far  from  cutoff,  3  3  (n  kQ-kx)  *  nkQ, 

thus  from  (4.2a)  and  (4.2b)  we  have  iEox/Eosi  =  3/kx  >>  1,  thus  we  will 

neglect  EQS  as  a  field  source  and  consider  our  straight  guide  modes  to  have 

only  x-directed  components.  Recalling  (2.49),  where  EQr  *►  Eqw  (i.e., 

E  becomes  E  in  our  local  coordinate  system  for  the  straight  guide  , 
or  ow 

E  ->  E  ,  and  E  ->  E  ,  we  conclude  that  for  the  special  case  we  are 
oq  os  oz  ox 

now  addressing  <t>m(x,w)  =  0  or,  in  particular 


*s(x,w)  =  0 


(C.  7) 


Rewriting  (1.4)  in  terms  of  the  local  coordinates  of  the  straight 
guide. 


o 

=  r 


tE„x6H-6ExH  } 

J  0  0 

-d 


a  dx 
w 


|x=b 


(C.8) 


C-3 


The  straight  guide  fields  and  H~  can  be  written  in  the  form 


Eo  =  £oxax 


(C.Sa) 


H  =  -  H  a  +  H  a 
o  ow  w  os  s 


(C.do) 


because  of  the  scalar  product  with  aw  that  occurs  in  (C.8)  and  the  forms 


of  E”  and  H“,  we  conclude  that  we  need  only  find  the  components  sE  of 


6E  and  <sHs  of  6H  in  order  to  evaluate  C. 


p 

Using  (C.7)  in  (1.5c)  and  recalling  that  we  have  defined  $  =  Ex 

and  <fm=  Hx  (where  the  subscript  z  becomes  x  and  the  subscript  0  becomes 
s  in  our  local  coordinate  system)  we  have 


ko"2H!  '  iWE  W"  ^  i  <— ?> 

3  2 


(C.10) 


where  we  have  inserted  the  Hankel  Transform  representation  of  Hf.  Substituting  from 

C  <> 

(1.10c)  in  C.10)  gives 


2  2s  •  3<{>< 

k0n2Ha  1We  Jf 


*1 4 


2  7  3  he 

J  — §  Ov(k  ar)da 
;0  0  3Z 


(C.n) 


However,  since  from  (3.23)  we  see  that  V*  is  linear  in  z,  it  follows 

32<S>®  324>® 

from  (2.49)  that  - f  -  0,  thus  — =  0,  and  (C. 11)  becomes 

32  3z 


hs  -  -r! 

4  ;0k0  3W  "s 


“  «c(X»W) 


(C.12) 


3  3  3 

where  we  have  used  —  =  ~ 

o*  Qj  d" 


If  we  now  substitute  (4.1)  into  (3.23)  we  see  that 


=  ^~k  0^(x,w) 
3w  2  o  s 


so  use  of  (C.13)  in  (C.12)  yields 


H*  A  6H.  =  «eU,w) 

O 


(C.13) 


(C. 14) 


With  (C. 9a),  (C.9b),  and  (C.14)  we  find  the  following: 


<■-  -  _ 2  0  (x,w)  E  (x,w) 

E  xoH'a  =  — —  s  ox 

owC 

o 


(C. 15) 


6EXH  'a  =  -  <b  (x,  w)  H  (x,w) 
o  w  s  o4 


Substitution  of  (C.15)  and  (C.16)  in  (C.8)  yields 


■  [(?) 


E  (x,w)+H  .  (x, w)>  dx 
ox  o4 


(C.16) 


(C.17) 


Finally,  we  can  use  (4.2b)  and  (4. 2d)  in  (C.17)  to  generate  the 


expression 


C  *  *  |  4  (x,b) 

J  s 

-d 


2.2  J2  i 
n  k  -k 


1  n  K  -K  K 

r  -S-*  “s  V™’ 

1  X 

k 

—  sin  k  (bin)  sin  k  (x+£)  dx 

k  w  I  x 

x  J  J 


n  k  k 


(C. 18) 


-•s*. 


C . 2 . 2  An  Explicit  Representation  of  0  _ (x,w) 


Having  (C.18),  our  objective  becomes  that  of  finding  an  explicit  form 
of  4's(x,w>  for  the  special  case  under  consideration.  We  begin  by 
rewriting  (3.23)  in  the  abbreviated  form  (in  terms  of  local  coordinates) 


VC  (w,w‘ ,x,x' )  =  C  Vs (w,w* ) u& (x,x* ) 
s  o  e  e 


(C. 19a) 


where 


c  t  2  _ !_ 

o  2  2,  „  A  k  R 

16tt  n  k  R  l  2  o  i 

loo  \ 


1  \3/2 

2  I  e 


(C. 19b) 


k  A  (w~b/2)  k  A  (w* -b/2) 
,.S  ,  ...  o  2  o  2 

V  (w ,  w  )  A  e  e 

e  = 


vC. 19c) 


u  (x, x 1 )  A  (k  x'-k)  (k  x-k) 
e  =  o  =  o  = 


(C. 19d) 


where 


2  2  2. 1/2 
r‘2(Vni) 


(C. 19e) 


Upon  substitution  of  (C.19a)  and  (4.2b)  into  (2.64)  we  find 

2  2  2 
2  n  k  -k 

4>^(w,x)  =  -AC  w (n  -n  )R  (k,+  2—r)  — r~-  — 

s  o  2  o  2  ax2  0kx 


)  V  (w,w*)  coz  k  (w'+  n)dw*  J  u  (x,x*)  sin  k  (x'+Odx*  (C.20) 

*  G  W  '  „  6  X 


The  integral  on  (0,b)  in  (C.20)  becomes,  upon  inserting  (C.19c), 


VS(w,w')  cos  k  (w'  +  n)dw '  = 
e  w 

k  A.(w-|)  fb  k  A  (w'  -b/2) 

e  °  J  e  °  cos  k  (w'  +  r^dw' 

*  w 

o 


(C. 21) 


Carrying  out  the  integration  on  the  RHS  of  (C.21)  and  using  (D-4a) 


(D-4b) »  (D-6a) ,  and  (D-6b)  yields 
b 


C.  A 

1  = 


k  X  (w'-b/2) 
o  2 


cos  k  (w'+n)dw’  = 
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^  x  k  fsgn[cos(kwb>]  sgn[k“-y2] 

koY2(n  " 1 1  2  t  L 


(k  bX  )/2  -(k  bX  )/2  -i 

o  2  o  2  I 


[‘ 


+  Y4  j sgn[cos(kwb) ] 


?  ,  (k  bX,)/2  -(k  bAJ/2 

f.  ^  o  2  o  2 

sgn  [k  -y. J  e  +  e 

w  4 


(C. 22) 


Combining  (C.22)  with  (C.21)  we  have 


V^(w,w‘)  cos  k  (w'+n)dw  =  C,  e 
e  w  1 


koX^  (w---) 


(C. 23) 


The  integral  on  [~d,0]  in  (C.20)  becomes,  upon  inserting  (C.19d), 


f  s 

u  (w,w' )  sin  k  (x'+4-)dx'  - 

J  e  x 


-d 


| 

(k  x-k)  J  (k  x'-k)  sin  k  (x'+F; dx' 
°  =  -d  °  =  x 


(C. 24 ) 


Evaluating  tne  integral  on  the  RHS  of  (C.24)  and  using  (D-7a) ,  (D-7b) , 
(4.6a),  and  (4.6b)  we  find 


r  r  * 

C  A  (k  x'-k)  sin  k  (x'+£;dx'  =  {  — 

2  =  J  o  x  '  k 

-d  '• 


(n^d+n2) 
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(n2kx+n  Y2 


iiiii  ii-Utl  iWii 
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C-7 


_5 _ i  ^  r 

>nS2)^2-l  +Rx  L( 


,42  4,  2 

<n  Y3+n2Kx) 


,  4  2  421/2 

(n  y  +n  k  ) 

1 3  lx 


+  sga  icos(k  d)]  sgn  l^n^^n4Y  Y  I  TT^TT/.'  .  (C*2:>) 

(n2kx+n  V  j 


Combining  (C.25)  with  (C.24)  we  have 


u  (x,x‘)  sin  k  (x'+C)dx’  =  C_  (k  x-k) 

J  e  x  2  o  = 

-d 

We  can  now  employ  (C.19b),  (C.23),  and  (C. 26)  with  (C.20)  to 


(C. 26) 


produce  the  result 


■gx)  =[  AC  C 

i. 


2  2  2  2  2  2 
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16tt  3k 
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k  X„(w~~)  k  R  q  (0) 
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e  e  (k 


x-k)  1 


(C. 27) 


where  we  have  also  used  we  =  - —  . 

°  ‘•o 


C.2.3  An  Explicit  Representation  of  C 


Having  (C.27),  we  need  only  use  it  with  (C.18)  to  obtain 


C  =  I  -Ac.^  sgn  [cosfk^b))  sgn  [k“— Y4 1 

k  k  (n2k2-k2) (n2-n2)1/2  fx  (n2k2-k2)  +  n2k  Y,1 
w  o  1  .  * _ 2 _  \_  2  ox  o  4J 


2  2  2 

16t23  kxX2C0 


k  R  X_ 
o  o  2 


<koi2b,/2  Voqi<0) 

e  c 


(C.23) 


D-l 


APPENDIX  D 

THE  EXPLICIT  REPRESENTATION  OF  TRIGONOMETRIC 
FUNCTIONS  IN  TERMS  OF  PROBLEM  PARAMETERS 

During  the  course  of  the  analysis  in  Appendix  C  we  employ  a  set  c-0 
identities  which  express  trigonometric  functions  explicitly  in  terms  of 
the  eigenvalues  K  ,  «w  and  the  material  parameters  of  our  problem.  Tne 
development  of  these  identities  from  equations  (4.4a),  (4.4b),  (4.6a), 
and  (4.6b)  is  a  straightforward  arithmetical  procedure,  but  one  must 
exercise  care  in  accounting  for  a  multitude  of  possible  choices  of 
numerical  sign  if  ambiguous  results  are  to  be  avoided. 

Consideration  of  (4.4a)  and  (4.4b)  reveals  that  the  eigenvalues 

K  ,  Kw  may  be  either  less  than  or  greater  than  zero.  Indeed,  review  of 

the  field  expressions  (4.2a)  through  (4.2f)  shows  that  the  numerical  signs 

of  Kx  and  i<w  don't  effect  these  formulae  at  all;  i.e.,  the  field 

expressions  depend  on  |Kx|  and  ,KW,.  We  are  thus  free  to  choose  the 

numerical  signs  of  K  and  K  and  will  choose  these  to  be  the  positive 

X  w 

solutions  of  the  eigenvalue  equations  (4.4a)  and  (4.b). 

Consider  now  equations  (4.6a)  and  (4.6b).  Since  we  choose  X  >  0  and 
Kw  >  0,  and  since  y^  >  0  and  y^  >  0,  we  can  always  find  solutions  £ 
and  n  such  that 


“77 

2K. 


<  £  <  0 


(D.la) 


D-2 


^ ' n " 0 


(D.lb) 


Choosing  Z  and  r,  as  in  (D.la)  and  (D.lb),  respectively  we  have 


sin(Kwn)  <  0 
cos(Kwn)  >  0 
sin(KxO  <  0 
cos(Kx^)  >  0 


(C.2a) 

(D.2b) 

(D.2c) 

(D.2d) 


We  can  now  proceed  to  outline  the  derivation  of  the  identities  we 
require  in  Appendix  C. 

Using  the  trigonometric  identites 


sin  6  = 


(1  +  tan26)i5 


sgnTcos  8  ] 
(1  +  tan2e) 


(0.3a) 


(D.3b) 


with  (4.6b)  and  (D.2a),  (D.2b),  we  find 


sin(K n)  - - -s - jn 

k0{"2  -  n2) 


(3.4a) 


cos(k  r,)  =  - o - ?r~r 

k0(n2  -  n2)** 


(D.4b) 


Similarly,  use  of  (4.4b)  with  (D.3a),  (D.3b)  gives 

sgn[cos(Kwb)]sgn[K^  -  ■y|3(2Kwy4) 
Sfn<K«b)  "  k2(n2  -  „2) 


(D.Sa) 


D-3 


cos(K  b) 
w 


sgnCcos(Kwb)] 1 -  y*\ 
-  n2) 


(0.5b) 


The  use  of  both  (D.4a),  (0.4b)  and  (0.5a),  (0.5b)  with  the  standard 
trigonometric  identities  for  sines  and  cosines  of  angular  sums  yields 


sin[Kw(b  +  n)]  =  sgn[cos(Kwb)]sgn[K2  -  y^J  - 2~2\ 

kQ(n  -  n2)J 


(D .6a) 


,2  2- 


cosLKw(b4n)]  -  sgn[cos(Kwb)]sgn[Kw  -y43  - \ 

k0(n  -  n2) '' 

Now,  proceeding  in  a  similar  way,  one  employs  (4.6a)  with  (0.3a), 
(0.3b)  and  (0.2c),  (D.2d)  to  produce 


(D.6b) 


Sln  ^  =  /A2  +  A2!15 

(n  y3  +  n^) 


C.7a) 


C°S  ^  =  / n4  2  ,  4,2} 

(n  y3  +  ^Kx) 


(D.7b) 


Using  (0.3a)  and  (0.3b)  with  (4.4a)  one  finds 


sin(Kxd)  = 


sgn[cos(Kxd)]sgn[n2]n2K2  -  n4Y2Y3]n2Kx(n2Y2  +  n2Y3) 


(n^K2  +  n4Y2)is(n4K2  +  n^2)^ 


(0.8a) 


sgn[cos(Kxd)]|n2n2K2  -  n4Y2Y3| 


COS(K  d)  =  — v— 3 - - - - TTir 

x  (n^  +  „\2)  („^  t 


(0.8b) 


Finally,  use  of  (0.7a),  (0.7b)  and  (0.8a),  (0.8b)  with  the  trigonometric 
formulae  for  angular  sums  and  differences  yields 


sin[K  (C  -  d)]  = 


2  2  2  4  2 

-  sgn[cos(Kxd)]sgn[n1n2Kx  -  n  y2y3]"2Kx 
(n2Kx  +  n^Z)h 


(D.9a) 


If 


cos[Kx(£-d)]  = 


sgn[cos(Kxd)]  sgn  [n^Kjj  -  n4Y2y3]n2y2 
(n2Kx  +  n4y2^ 


(D.9b) 


\ 


cos  Kw(w+rydw 


Rewriting  (1.4)  in  terms  of  the  local  coordinates  or  tne  siraiym. 


guide. 


{E~x5H-6ExH~ } 


awdx 


x=b 


(C.8) 


